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ABSTRACT 

We use new multi- wavelength radio observations, made with the VLA and Effelsberg 
telescopes, to study the magnetic field of the nearby galaxy M51 on scales from 200 pc to sev- 
eral kpc. Interferometric and single dish data are combined to obtain new maps at AA3, 6 cm 
in total and polarized emission, and earlier A20 cm data are re-reduced. We compare the spa- 
tial distribution of the radio emission with observations of the neutral gas, derive radio spectral 
index and Faraday depolarization maps, and model the large-scale variation in Faraday rota- 
tion in order to deduce the structure of the regular magnetic field. We find that the A20 cm 
emission from the disc is severely depolarized and that a dominating fraction of the observed 
polarized emission at A6 cm must be due to anisotropic small-scale magnetic fields. Taking 
this into account, we derive two components for the regular magnetic field in this galaxy: the 
disc is dominated by a combination of azimuthal modes, m = + 2, but in the halo only an 
m = 1 mode is required to fit the observations. We disuss how the observed arm-interarm 
contrast in radio intensities can be reconciled with evidence for strong gas compression in the 
spiral shocks. In the inner spiral arms, the strong arm-interarm contrasts in total and polar- 
ized radio emission are roughly consistent with expectations from shock compression of the 
regular and turbulent components of the magnetic field. However, the average arm-interam 
contrast, representative of the radii r > 2 kpc where the spiral arms are broader, is not com- 
patible with straightforward compression: lower arm-interarm contrasts than expected may 
be due to resolution effects and decompression of the magnetic field as it leaves the arms. We 
suggest a simple method to estimate the turbulent scale in the magneto-ionic medium from the 
dependence of the standard deviation of the observed Faraday rotation measure on resolution. 
We thus obtain an estimate of 50 pc for the size of the turbulent eddies. 

Key words: galaxies: spiral - galaxies: magnetic fields - galaxies: ISM - galaxies: individual: 
M51 



1 INTRODUCTION 

The Whirlpool galaxy, M51 or NGC5194, is one of the classical 
grand-design spiral galaxies. The two spiral arms of M51 can be 
traced through more than 360° in azimuthal angle in numerous 
wavebands. M51 is probably perturbed by a recent encounter with 
its companion galaxy NGC 5195. Such interactions usually result in 
enhanced star formation, either localized or global, as tidal forces 
and density waves compress the interstellar medium. In M51 this 
may have resulted in two systems of density waves { [Elmegreen et| 
|al.|1989| ). 

M5 1 was the first external galaxy where polarized radio emis- 
sion was detected ( [Mathewson et al.|1972| ) and one of the few ex- 
ternal galaxies where optical polarization has been studied ( Scarrott 
let al.|1987| ). rNeininger|(T992l ) and |Horellou etaL| ( [T992| found that 
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the lines of the regular magnetic field in M51 have a spiral shape, 
but their resolution was too low to determine how well the field is 
aligned with the optical spiral arms. |Horellou et al. ( 1992 ) realized 
that, at wavelengths of A ^ 18 cm, only polarized emission from 
a foreground layer reaches the observer because of Faraday depo- 
larization, and the Faraday rotation measures obtained at the longer 
wavelengths are much smaller than those at shorter wavelengths. 
[Heald et ar| ( |2009| ) observed a fractional polarization at A22 cm of 
around 5% in the optical disc, increasing to 30% at large radii, and 
Faraday rotation measures of a similar magnitude to those found 
by [Horellou et"ar] ( |1992| ). [Berkhuijsen et al.| \1991) analyzed the 
global magnetic field of M51 using radio data at four frequencies 
and found that the average orientation of the fitted magnetic field 
is similar to the average pitch angle of the optical spiral arms mea- 
sured by [Howard & Byrd| (19901 ). [Berkhuij sen et al.| ( [T997t repre- 
sented the magnetic field in the disc of M51 as a superposition of 
periodic azimuthal modes, with about equal contribution from the 
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axisymmetric m = and the bisymmetric m = 1 ones. Their fit 
contains a magnetic field reversal at about 5 kpc radius which ex- 
tends over a few kpc in azimuth. Furthermore, Berkhuijse n et al.] 
( 1997 ) found evidence for an axisymmetric magnetic field in the 
halo of M51 (visible at A18 cm and A20 cm) with a reversed direc- 
tion (inwards) with respect to the axisymmetric mode of the disc 
field (outwards). 



Table 1. The merged VLA and Effelsberg maps discussed in this paper, 
their half-power beam widths, imaging weighting schemes and their r.m.s. 
noises in total and polarized emission, aj and a pj , respectively 



The results of| Berkhuijsen et ah] ( | 1997) ) indicate that the mag- 
netic field of M51 is strong and partly regular, with some inter- 
esting properties. However, some of their results may be affected 
by the low angular resolution which was limited by the Effelsberg 
single-dish data at A2.8 cm. Furthermore, no correction for missing 
large-scale structure could be applied to the VLA A6 cm polariza- 
tion data because no single-dish map was available at that time. 



Here we present a refined analysis based on new data with 
higher resolution and better sensitivity. Our new surveys of M51 
in total and linearly polarized A3 cm and A6 cm radio continuum 
emission combines the resolution of the VLA with the sensitivity 
of the Effelsberg single-dish telescope. The new maps are of com- 
parable resolution to the maps of the CO(l-O) emission ( Heifer et 
aTpOOS j and mid-infrared dust emission ( S auvage et al. 1996 . Re- 



gan et al.|2006| ) and are of unprecedented sensitivity. The shape of 



the spiral radio arms and the comparison to the arms seen with dif- 
ferent tracers were discussed in a separate paper ( ^Patrikeev et al.| 
[2006] ). 



The spiral shocks in M51 are strong and regular ( [Aalto et aT] 
|1999j and offer the possibility to compare arm-interarm contrasts 
of gas and the magnetic field. We analyse our new data to try to 
separate the contribution to the observed polarized emission from 
regular (or mean) magnetic fields and anisotropic random magnetic 
fields produced by compression and/or shear in the spiral arms. 



The new high-resolution polarization maps allow us for the 
first time to investigate in detail the interaction between the mag- 
netic fields and the shock fronts. Results from the barred galax- 
ies NGC 1097 and NGC 1365 showed that the small-scale and 
the large-scale magnetic field components behave differently, i.e. 
the small-scale field is compressed significantly in the bar's shock, 
while the large-scale field is hardly compressed ( Beck et al.|2005| ). 
This was interpreted as a strong indication that the large-scale field 
is coupled to the warm, diffuse gas which is only weakly com- 
pressed. We investigate whether a similar decoupling of the reg- 
ular magnetic field from the dense gas clouds is suggested by the 
observed arm-interarm contrast in M5 1 . 
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2 OBSERVATIONS AND DATA REDUCTION 
2.1 VLA observations 

M51 was observed in October 2001 at A3. 5 cm with the VLaE] m 
the compact D configuration and in August 2001 at A6.2 cm in the 
C configuration. Two pointings of the array, at the northern and 
southern parts of the galaxy, were required to obtain complete cov- 
erage of M51 at A3. 5 cm. The data were edited, calibrated and im- 
aged using standard AIPS procedures and VLA calibration sources. 
After initial examination of the data and the flagging of bad visibil- 
ities, self-calibration was used. After correction for the pattern of 
the primary beam the two pointings at A3. 5 cm were mosaiced us- 
ing the AIPS task LTESS. The calibrated A6.2 cm C-array uv-data 
were combined with existing A6.2cm D-array data ( [Neininger &| 
|Horellou|1996l ). 

A new A20 cm map in total power and polarization is pre- 
sented in Fig. [3] based on the C-array data of , Neininger & Horel-| 
[Iou| ( |1996| ) combined with D-array data of |Horellou et al.| ( |1992| ) 
All of the observations were re-reduced for this work, the two data 
sets were combined and then smoothed to a resolution of 15^^ to 
improve the signal-to-noise ratio. The resulting maps at A20 cm 
contain information on all scales down to the beamsize as the pri- 
mary beam of the VLA at A20 cm in the D-array is ^ 30^ about 
twice the size of M5 1 on the sky. 

Different weighting schemes were used in the final imaging 
of the data to produce maps with either high resolution (natural 
weighting) or high signal-to-noise (uniform weighting) and for a 
compromise between the two extremes (robust weighting, obtained 
by setting the parameter ROBUST=0 in the AIPS task IMAGR). 
Maps of the Stokes parameters /, Q and U were produced in each 
case. The optimum number of iterations of the CLEAN algorithm 
used to produce the images was determined individually for each 
map. Following slight smoothing, the Q and U maps were com- 
bined to give the polarized intensity PI = y^Q^ + [7^, using 
a first-order correction for the positive bias ( (Wardle & Kronberg] 
17974) ). 



The basic parameters we adopt for M51 are: centre's right 
ascension 13^29" ^52!709 and declinati on (J2000) +47° 11^4 2^59 
( |Ford et al.|l'985l ); distance 7.6 Mpc, (Ciar dullo et al.|2002 ), thus 
1" ^ 37 pc; position angle of major axis —10° (0° is North) 
and inclination —20° (0° is face-on) ( |Tully||1974| ). The inclina- 
tion is measured from the galaxy's rotation axis to the line-of- sight, 
viewed from the northern end of the major axis, and its sign is im- 
portant for the geometry of the model discussed in Sect. [6] 



2.2 Effelsberg observations and merging 

In order to correct the VLA maps at A3 cm for missing extended 
emission we made observations of M51 in total intensity and po- 



1 The VLA is operated by the NRAO. The NRAO is a facility of the Na- 
tional Science Foundation operated under cooperative agreement by Asso- 
ciated Universities, Inc. 
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Figure 1. (a) A3 cm and (b) A6 cm radio emission at Ih" resolution from VLA and Effelsberg observations, overlaid on a Hubble Space Telescope optical 
image (image credit: NASA, ESA, S. Beckwith (STScI), and The Hubble Heritage Team (STScI/AURA)). Total intensity contours in both maps are at 
6, 12, 24, 36, 48, 96, 192 times the noise levels of 20 /iJy/beam at A3 cm and 30 /xJy/beam at A6 cm. (Note that the roughly horizontal contours at the left 
edge of panel (a) are artefacts arising from mosaicing the two VLA pointings.) Also shown are the S-vectors of polarized emission: the plane of polarization 
of the observed electric field rotated by 90° , not corrected for Faraday rotation, with a length proportional to the polarized intensity {PI) and only plotted 
where PI ^ 3crp/. 



larization with the 100-m Effelsberg telescope [^in December 2001 
and April 2002 using the sensitive A3. 6 cm (1.1 GHz bandwidth) 
receiver. We obtained 44 maps of a 12^ by 12^ field around M51, 
scanned in orthogonal directions. Each map in /, Q and U was 
edited and baseline corrected individually, then all maps in each 
Stokes parameter were combined using a basket weaving method 
( [Emerson & Grave, 1988 J . The r.m.s. noise in the final maps, after 
slight smoothing to 90^^ is 200/iJy/beam in total intensity and 
20 /iJy/beam in polarized intensity. 

At A6.2 cm, ten maps of a 41^ by 34^ field were observed in 
November 2003 with the 4.85 GHz (500 MHz bandwidth) dual- 
horn receiver. The combined data resulted in a new 180^^ resolution 
image with r.m.s. noise of 250/iJy/beam in total intensity and 
25 /iJy/beam in polarization. 

Maps from the VLA and Effelsberg were combined using the 
AIPS task IMERG. A useful description of the principles of merg- 
ing single dish and interferometric data is given by |Stanimirovic| 
( |2002| ). The range of overlap in the uv plane between the two im- 
ages (parameter uvrange in IMERG) was estimated as follows: 
we assumed an effective Effelsberg diameter of about 60 m to esti- 
mate the maximum extent of the single dish in the uv space (1.7 kA 
at A3. 6 cm and 1.0 k A at A6.2cm) and used the minimum sep- 
aration of the VLA antennas in the D-array configuration, 35 m, 



^ The Effelsberg 100-m telescope is operated by the Max-Planck-Institut 
fiir Radioastronomie on behalf of the Max-Planck-Gesellschaft. 



to calculate the minimum coverage of the interferometer in the uv 
space (1.0 k A at A3. 6 cm and 0.6 k A at A6.2cm). We then varied 
these parameters in order to find the optimum overlap in the uv 
space by comparing the integrated total flux in the merged maps 
with that of the single dish maps; the optimal uv-ranges for merg- 
ing were found to be 1.0 ^ 1.6 kA at A3. 6 cm and 0.5 0.7 kA 
at A6.2 cm. Merging of the Q and U maps was carried out using 
the same optimum uv range as found for /. 

The fraction of total emission (Stokes /) present in the VLA 
maps - those produced using natural weighting, and hence with the 
highest signal-to-noise ratio - is about 30% at A3 cm and close to 
50% at A6 cm compared to the merged maps. Small-scale fluctu- 
ations in Q and U due to variations in the magnetic field orien- 
tation and Faraday rotation in M5 1 mean that the polarized emis- 
sion is less severely affected by missing large-scales (alternatively, 
the single dish detects the large-scale emission missed by an in- 
terferometer but simultaneously suffers from stronger wavelength- 
independent beam depolarization). At A3 cm the VLA map con- 
tains about 75% of the polarized emission present in the merged 
map, with about 85% present at A6 cm. 

Following the merging, the maps in /, Q and U were con- 
volved with a Gaussian beam to give a slightly coarser resolution 
and higher signal-to-noise ratio. The maps that are discussed in this 
paper are listed in Table [T] along with the r.m.s. noises in total and 
polarized intensity. 
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Figure 2. (a) A3 cm and (b) A6 cm (right) polarized radio emission at IW resolution from VLA and Effelsberg observations, overlaid on the same optical 
image as in Fig.[^ Polarized intensity contours in both maps are at 3, 9, 15, 21 times the noise level of 8 /iJy/beam at A3 cm and 10 /iJy/beam at A6 cm. 
Also shown are the B-vectors of polarized emission: the position angle of the polarized electric field rotated by 90°, not corrected for Faraday rotation, with 
the length proportional to the polarized intensity PI and only plotted where P ^ Sap/. 



3 THE M51 MAPS 

3.1 The spiral arms of M51 as seen in radio emission 

Figure [T] shows the total radio continuum emission and the B- 
vectors of polarized emission (the observed plane of linear polar- 
ization rotated by 90°) at A3 cm and A6cm, overlaid on a Hub- 
ble Space Telescope optical image. At the assumed distance of 
7.6 Mpc, the 1^" resolution corresponds to 560 pc and 590 pc 
along the major and minor axes, respectively. The distribution of 
polarized emission at A A3, 6 cm is shown in Fig.|2]at 15^^ resolu- 
tion. The extensive A20 cm total emission disc is shown in Fig. [3] 
also at 15^^ resolution. 

The total emission at A3 cm and A6 cm in Fig.[T]shows a close 
correspondence with the optical spiral arms whereas the A20 cm 
total emission (Fig. [3]) and A3 cm and A6 cm polarized emission 
(Fig.[2| are spread more evenly across the galactic disc. Compact, 
bright peaks of total emission coincide with complexes of H II re- 
gions in the spiral arms, as expected if thermal bremsstrahlung is a 
significant component of the radio signal at these peaks. The flatter 
spectral index {a < 0.6 where / oc z/~". Fig.]?]) and the absence of 
the corresponding peaks in polarized radio emission (Fig. |2]) sug- 
gest that a significant proportion of the centimeter wavelength ra- 
dio emission in these peaks is thermal. The extended A20 cm total 
emission and A3 cm and A6 cm polarized emission accurately trace 
the synchrotron component of the radio continuum. 

Ridges of enhanced polarized emission are prominent in the 
inner galaxy; some are located on the optical spiral arms but others 
are located in-between the arms. A detailed analysis of the location 



and pitch angles of the spiral arms traced by different observations 
and a comparison of the pitch angles with the orientation of the 
regular magnetic field is given by Patrikeev et al.|(|2006|). They find 
systematic shifts between the spiral ridges seen in polarized and to- 
tal radio emission, integrated CO line emission and infrared emis- 
sion, which are consistent with the following sequence in a den- 
sity wave picture: firstly, shock compresses gas and magnetic fields 
(traced by polarized radio emission), then molecules are formed 
(traced by CO) and finally thermal emission is generated (traced by 
infrared) . [Patrikeev et al . | ( [2006| also show that while the pitch an- 
gle of the regular magnetic field is fairly close to that of the gaseous 
spiral arms at the location of the arms, the magnetic field pitch an- 
gle changes from the above by around ±15° in interarm regions. 



3.2 The connection between polarized radio emission and 
gaseous spiral arms 

In Fig. |4] the A6 cm polarized emission and the orientation of the 
regular magnetic field in the central 8 kpc of M51 are overlaid onto 
an image of the spiral arms as traced by the CO(l-O) integrated line 
emission. Part of the polarized emission appears to be concentrated 
in elongated arm-like structures that sometimes coincide with the 
gas spiral. 

The correspondence between polarized and CO arms is good 
along most of the northern arm in Fig. |4] (called Arm 1 in the rest 
of the paper) and in the inner part of the southern arm (called 
Arm 2). Arm 1 continues towards the south where the polarized 
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Figure 3. Contours of A20 cm total radio emission at 1^" resolution, over- 
laid on the same optical image as in Fig.[T] Total intensity contours are at 6, 
12, 24, 36, 48, 96, 192 times the noise level of 20 /iJy/beam. Also shown 
are the S-vectors of polarized emission: the plane of polarization of the ob- 
served electric field rotated by 90° , not corrected for Faraday rotation, with 
a length proportional to the polarized intensity PI and only plotted where 
PI ^ 3crpj. 



emission longer coincides with the CO and optical arms, but be- 
comes broader further out (Fig.[2|. 

Moving along Arm 2, the excellent overlap of the radio po- 
larization and CO in the inner arm ends abruptly further in the 
south (Fig. |4]), beyond which the distribution of polarized emis- 
sion is very broad and located on the galactic centre side of the 
gaseous spiral arm: the peak of this interarm polarized emission 
corresponds very closely with the pronounced peaks in radial ve- 
locity, Vr — lOOkms"^, derived by Shetty et al. (2007 ) (shown 
in their Fig. 14). This may be an indication that strong shear in 
the interarm gas flow is producing this magnetic feature. At about 
6 kpc radius, the polarization Arm 2 crosses the CO and optical 
arms in the east, followed in the northeast by a bending away from 
the optical arm towards the companion galaxy NGC5195 (Fig. [2]). 
The northern part of Arm 2 is weakly polarized and rather irregular 
in total and CO emission. The whole space between the northern 
Arm 2 and the companion galaxy is filled with highly polarized 
radio emission (typically 15% at A6 cm). Arm 2 becomes well or- 
ganized again at larger radii (located at the western edge of Fig.[2|, 
where the total radio, polarized radio and CO emission perfectly 
coincide. 

West of the central region, between Arms 1 and 2 in Fig. [4] 
another polarization feature emerges which appears similar to the 
magnetic arms observed e.g. in NGC 6946 ( B eck & Hoernes|1996| ). 
However, in contrast to NGC 6946, Faraday rotation is not en- 
hanced in the interarm feature of M51 (see Fig. [9]). Some peaks of 
polarized emission between Arms 1 and 2 in the south and south- 
east (see a low-resolution image of Fig. [2]) and may indicate the 



Figure 4. Contours of the A6 cm polarized radio emission (VLA and Ef- 
felsberg combined) in the central ~ 3' x 4' of M5 1 at 8^^ resolution, along 
with Faraday-rotation corrected B -vectors, overlaid on the map of integrated 
CO(l-O) line emission of HelfCTeF^ (2003). Contours are at 3, 5, 7, 9 
times the noise level of 10 /iJy/beam. 



outer extension of this magnetic arm. Inside of the inner corotation 
radius, located at 4.8 kpc ( Elmegreen et al.|1989 ), this phenomenon 
can be explained by enhanced dynamo action in the interarm re- 
gions ( |Moss|1998||Shukurov|1998l|Rohde et al.|1999| ). 



3.3 Polarized radio emission from the inner arms and central 
region 

In the CO and Ha line emissions (Fig. |4] and the red regions in 
Fig. |5]), the spiral arms continue towards the galaxy center. The 
high-resolution CO map by Aalto et aL] (|1999) shows that the 
arms are sharpest and brightest between about 2^" and 50^^ dis- 
tance from the center. The arms become significantly broader and 
less pronounced inside a radius of about 0.8 kpc: this is inside 
the inner Lindblad resonance of the inner density-wave system at 
r ^ 1.3 kpc identified by Elm egreen et al.| ( p"989| t. 

The polarized emission at A" resolution (see Fig [6} is also 
strongest along the inner arms 1-2 kpc distance from the center, 
with typically 20% polarization. The arm-interarm contrast is at 
least 4 in polarized intensity (this is a lower limit as the interarm 
polarized emission is below the noise level at this resolution and we 
take a PI as an upper limit for the interarm value), larger than that 
of the outer arms, and is consistent with the expectations from com- 
pression of the magnetic field in the density-wave shock (Sect. [7]). 
The contrast weakens significantly for r < 0.8 kpc. This may be 
an indication that the inner Lindblad resonance of the inner spi- 
ral density wave is at r ^ 0.8 kpc rather than r ^ 1.3 kpc (as 
located by [Elmegreen et aL] ( |1989| )): the shock is probably weak 
around the inner Lindblad resonance. In total intensity, the typical 
arm-interarm contrast for the region of the inner arms is about 5. 
The actual contrast in the M51 disc alone may be stronger than this 
if there is significant diffuse emission in the central region from a 
radio halo, but this effect is hard to estimate. 

In the central region, two new features appear in polarized in- 
tensity which are the brightest in the entire galaxy (Fig. [6]). The 
first is a region 11^^ north of the nucleus with a mean fractional po- 
larization of 10% and an almost constant polarization angle. This 
feature coincides with the ring-like radio cloud observed in total 
intensity at A6 cm and at 1" resolution by Ford et al. ( |1985| ) who 
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Figure 5. (a) A3 cm and (b) A6 cm total radio emission and B-vectors (VLA and Effelsberg combined) in the central 3^ x 4' of M51 at 8^^ resolution overlaid 
on a Hubble Space Telescope image (http://heritage.stsci.edu/2001/10, image credit: NASA and The Hubble Heritage Team (STScI/AURA)). Contours are at 
6, 12, 24, 36, 48, 96, 192 times the noise level. B-vectors, not corrected for Faraday rotation, are plotted where PI ^ 3crp/. 
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Figure 6. A6 cm total radio emission (VLA and Effelsberg combined) in 
the central 1.4^ x 1.4^ of M51 at 4'' resolution. The B-vectors (not cor- 
rected for Faraday rotation) are shown where PI ^ ^crpi- 



also detected polarization in this region. The polarized emission in- 
dicates that the plasma cloud expands against an external medium 
and compresses the gas and magnetic field. The second feature of 
similar intensity in polarization is a ridge located along the east- 
em edge of the first region, extending east of the nuclear source, 
with 15% mean polarization and a magnetic field almost perfectly 
aligned along the ridge. Field compression is apparently also strong 
in this ridge. The nuclear source itself appears unpolarized at this 
resolution. 



4 SPECTRAL INDEX AND MAGNETIC FIELD 
STRENGTH 

Spectral index maps of the total radio emission at 15^^ resolution are 
shown in Fig.[7] The spectral indices are calculated using combined 
Effelsberg and VLA maps that contain signal on all scales down to 
the beamsize. There is therefore no missing flux in these maps due 
to the short-spacings problem of interferometers. 

Although there is good general agreement between the two 
spectral index maps, the spectral index between AA20, 3 cm is gen- 
erally slightly flatter than that between AA20, 6 cm; this is partic- 
ularly noticeable in the spiral arms. This is due to the lower res- 
olution of the Effelsberg A6cm map (180^^ beam against 90^^ at 
A3 cm): 180^^ is about the radius of the M51 disc. The coarse 
resolution made it more difficult to fine-tune the merging of the 
VLA and Effelsberg data in order to match the integrated fluxes 
in the single-dish and merged maps. This in turn led to a slight 
under-representation of the single-dish data at A6 cm in the merged 
map and thus to a slightly steeper spectral index. We believe the 
AA20, 3 cm spectral index, shown in Fig. [7^, to be more reliable. 

In both spectral index maps, one can clearly distinguish arms 
and interarm regions. The spectral index is typically in the range 
—0.9 ^ a ^ —0.6 in the spiral arms and — 1.2 ^ a ^ —0.9 in the 
interarm zones. Since the spectral index in the arms is not character- 
istic of thermal emission {a / —0.1) the arm emission must com- 
prise a mixture of thermal and synchrotron radiation. The flatter 
arm spectral index can then be explained by two factors: stronger 
thermal emission in the arms due to recent star formation and the 
production of H II regions and energy losses of cosmic ray elec- 
trons as they spread into the interarm from their acceleration sites 
in the arms. 

We can use the observed steepening of the spectral index of 
the total radio emission a to estimate the diffusion coefficient of 
cosmic ray electrons in M51. We assume that the sources of the 
electrons are supernovae in the arms, that the initial spectral index 
is ttsyn ^ — 0.5 and that the radio emission in the interarm region 
is predominantly synchrotron, so that asyn — —1.1 between the 
arms. A difference in the spectral indices of Aa = 0.5 is expected 
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Figure 7. The spectral index maps between (a) A20 cm and A3 cm and (b) A20 cm and A6 cm at Ih" resolution. Also shown are contours of total radio 
emission at A3 cm (left) and A6 cm (right); the contour lines are drawn at 6, 12, 24, 36, 48, 96, 192 times the noise levels of 20 /iJy/beam and 30 /iJy/beam 
at A3 cm and A6 cm, respectively. Spectral indices are only calculated where the signal at both wavelengths is 6 times the noise level. The error in the fitted 
spectral index is typically ±0.02 in the inner galaxy and spiral arms and ±0.06 in the interarm regions at radii > l^ The large red area at the left edge of (a) 
is due to increased noise in the overlap region between the two VLA-pointings at A3 cm. 



if the main mechanisms of energy losses for the electrons are syn- 
chrotron emission and inverse Compton scattering ( ^Longair|1994| ). 

In the next section we estimate magnetic field strengths of 
around 20 /iG for the inter- arm regions. In these magnetic fields, 
an electron emitting at 5 GHz has an energy of about 4 GeV. We 
can estimate the lifetime of cosmic ray electrons emitting the syn- 
chrotron radiation as 

^16 ^tot± 



( |Lang|p^99| Sect. 1.25) where the frequency viq is measured in 
units of 16 MHz, the total magnetic field strength in the plane of the 
sky Btot± is measured in /xG and we have taken Btot± ^15 /xG 
in the interarm regions. Taking L = 1 kpc as the typical distance a 
cosmic ray electron travels from its source in a supernova remnant 
to the inter-arm region, yields the diffusion coefficient D of the 
electrons, 
r2 



D 



4 X 10^^ cm^s"\ 



which is compatible with the value of D 1-10 x 10^^ cm^ 
estimated by Strong & Mosk alenko|(T998 ) for the Milky Way. Note 
that cosmic ray electrons producing radio emission at cm wave- 
lengths, propagating for a few kpc in /xG strength magnetic fields, 
give diffusion coefficients in this range: our estimate for D is not a 
unique property of the cosmic rays in M5L 



4.1 Total magnetic field 

In order to derive the strength and distribution of the total magnetic 
field, we have made a crude separation of the A6 cm map at 1^" 
into its nonthermal /syn and thermal hh components. We assumed 
that the thermal spectral index is everywhere ath = —0.1 as ex- 
pected at cm wavelengths (e.g. |Rohlfs & Wilson|[l999| and that 
the synchrotron spectral index is asyn = —1.1 everywhere, as ob- 
served in the interarm regions (Fig. [7^). We are constrained in our 
choice of asyn by two considerations: if agyn > —1.1 then ther- 
mal emission is absent from the whole inter-arm region, whereas 
H-a emission is detected; if ttsyn < —1.1 we find that the degree 
of polarization approaches is maximum theoretical value of 70% 
in many regions, which is implausible for our resolution of 570 pc. 
The average thermal emission fraction at A6 cm is 25%. 

Assuming equipartition between the energy densities of the 
magnetic field and cosmic rays, a proton-to-electron ratio of 100 
and a pathlength through the synchrotron-emitting regions of 1 kpc, 
estimates for the total field strength are shown in Fig. |8] applying 
the revised formulae by |Beck & Krause| ( |2005| ). 

The strongest total magnetic fields of about 30 /xG are ob- 
served in the central region of M51. The main spiral arms host to- 
tal fields of 20-25 /iG, while the interarm regions still reveal total 
fields of 15-20 /xG. This is significantly larger than in spiral galax- 
ies like NGC 6946 ( |Beck|2007) and M33 ( |Tabatabaei et al.|2008| ). 
These two galaxies have similar star-formation rates per unit area 
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Figure 8. Total magnetic field strength, derived from the A6 cm emission 
assuming equipartition between the energy densities of magnetic fields and 
cosmic rays (colour scale in /xG) along with contours of neutral gas den- 
sity (a combination of CO jHelfer et al.[2003) and HI jRots et al.|1990t 
observations) plotted at 1, 4, 8, 16, 32, 64% of the maximum value. 



as M51, but weaker density waves, so that compression is probably 
higher in M5 1 . 



4.2 Ordered magnetic field 

The strength of the ordered magnetic field can be estimated from 
that of the total field using the degree of polarization. This method 
gives field strengths of 11-13 /iG in the inner spiral arms, 8- 
10 /xG in the outer spiral arms and 10-12 /xG in the inter- arm re- 
gions. However, these values can only be attributed to a regular 
(or mean) magnetic field if the unresolved random component is 
purely isotropic (see'Sokoloff et al.ll998l Sect. 5.1). The observed 
maximum degree of polarization of around 40% can equally be pro- 
duced by an anisotropic random field whose degree of anisotropy 
is about 2: that is if the standard deviation of the fluctuations in 
one direction on the plane of the sky is twice as large as in the or- 
thogonal direction. In Sect. [6] we shall see that there is only a weak 
signature of a regular field in the observed multi-frequency polar- 
ization angles and that most of the polarized emission does indeed 
arise due to anisotropy in the random field. 

The difference between the total and ordered magnetic field 
strengths gives an estimate for the isotropic random magnetic field 
strength of 18 /xG, or 1.5 times the ordered field, in the arms and 
13 fiG, or 1.2 times the ordered field, in the inter-arms. If the main 
drivers of (isotropic) turbulence are supernova remnants, then the 
preferential clustering of Type II supernovae in the spiral arms is 
compatible with the higher fraction of isotropic random field in the 



arms. So a significant fraction of the magnetic field consists of a 
random component that is isotropic on scales less than 500 pc. 

4.3 Uncertainties 

Our assumption of a synchrotron spectral index that is constant 
across the galaxy is crude and an oversimplification, even though 
the value used, asyn = —11, can be somewhat constrained by 
other data as described above. We would expect that ttsyn should 
be closer to —0.5, the theoretical injection spectrum for electrons 
accelerated in supernova remnants, in parts of the spiral arms. This 
limitation results in an overestimate (underestimate) of the thermal 
(nonthermal) emission in the arms. 

In principle we could combine the data at all three frequencies 
and simultaneously recover /gyn, hh and asyn at each pixel. How- 
ever we defer a more robust calculation, interpretation and discus- 
sion to a later paper. 

The equipartition estimate depends on the input parameters 
with a power of only 1/(3 + ttsyn) ~ 0.24, so that even large 
input errors hardly affect the results. Further errors are induced by 
the underestimate of the synchrotron emission in the spiral arms by 
the standard separation method (see above). In M33, |Tabatabaei et| 
[aTlpOOVj ) found that the standard method undestimates the average 
nonthermal fraction by about 25%. In star-forming regions of the 
spiral arms, the nonthermal intensity can be a factor of two too 
small, which leads to an equipartition field strength which is 20% 
too low. In the same regions, the synchrotron spectral index is too 
steep by about 0.5 which overestimates the field strength by 15%. 
Interestingly, both effects almost cancel in M33. 

The equipartition assumption itself is subject of debate. 
Equipartition between cosmic rays and magnetic fields likely does 
not hold on small spatial scales (e.g. smaller than the diffusion 
length of cosmic rays) and on small time scales (e.g. smaller than 
the diffusion time of cosmic rays). 



5 FARADAY ROTATION AND DEPOLARIZATION 
5.1 Faraday rotation 

The nonthermal radio emission from the arms has a relatively low 
degree of polarization (typically 25% at A6 cm and 15^^ resolution) 
so that unresolved, tangled or turbulent, magnetic field dominates 
in the arms. In contrast, the interarm regions are up to 40% po- 
larized and host a significant fraction of magnetic fields with ori- 
entation ordered at large- scales. Whether these fields are coherent 
(regular) or incoherent (anisotropic turbulent) can be decided only 
with the help of Faraday rotation measures. 

One might expect from the well-ordered, large-scale spiral 
patterns of the polarization vectors of Figs. [2] and [5] that the reg- 
ular magnetic field would produce an obvious pattern in Fara- 
day rotation. (See the rotation measure map of M3 1 in Fig. 1 1 of 
[Berkhuijsen et al.| ( [20 03 ) for an example of clear rotation measure 
signal arising from a well-ordered magnetic field.) However, the 
Faraday rotation measure map shown in Fig. [9] is dominated by 
strong fluctuations in rotation measure, with a magnitude of order 
100 rad m~^. We apparently have a paradox: the orientation of the 
regular magnetic field follows a systematic spiral pattern on scales 
exceeding 1 kpc but it does not produce any obvious large scale 
pattern in RM. Even the magnetic arms located between the CO 
arms do not immediately exhibit any large-scale non-zero RM. If 
the ordered field seen in polarized intensity with an equipartition 
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Figure 9. Rotation measures between A A3, 6 cm, at 15" resolution, over- 
laid with contours of Ho; emission ( Greenawalt et al.|1998| at the same 
resolution, plotted at 4, 8, 16, 32% of the map maximum. Data were only 
used where the signal-to-noise ratio in polarized intensity exceeds three. 
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Figure 10. Distribution of rotation measures between A A3, 6 cm, shown in 
Fig.[9| note that the map is oversampled and so the histogram of pixel counts 
does not represent statistically independent data points. Data were only used 
where the signal-to-noise ratio in polarized intensity exceeds three. Solid 
line is the best-fitting Gaussian to the histogram. 



strength of about 10 /iG (Sect. [4]) was fully regular, we would have 
|RM| ^ 700 rad near the major axis of M51 with a system- 
atic decrease moving away from this axis in azimuth, which is not 
observed (we adopt (rie) = 0.1 cm~^, h = 400 pc and an incli- 
nation of 20° for this estimate, see Sect. \6.2\ . Note also that the 
overlaid Ha contours in Fig. |9] are not generally coincident with 
regions of strong rotation measure. 



Observational uncertainties in the measured Stokes parame- 
ters can be a source of the fluctuations in Fig. |9] The uncertainty 
in RM between A3 cm and A6 cm, denoted here ^RM, depends 
on the signal-to-noise ratio in the polarized intensity, S, and the 
difference in the squares of the wavelengths ^A^ as 



rr^^/r 1 280 ^ _2 

^RM=^ ^— — radm . 

V2S^A2 ^ 



(1) 



Due to the steep spectral index of the synchrotron emission and 
the weak Faraday depolarization between A3 cm and A6 cm, S is 
lower at A3 cm and we use these values to estimate ^RM. The RM 
fluctuations from the noise in the observed polarization signal are 
±10 rad m"2 (for typical S > 30 in the central r < 90'' of Fig.jo}. 
RM maps at 8'' resolution have lower H and are dominated by 
noise fluctuations, although some strong rotation measures, above 
the noise level, are also present. 

The structure function of the RM fluctuations is flat on scales 
up to 3', whereas RM fluctuations due to the Milky Way fore- 
ground in the direction of M5 1 in the model of | Sun & Reich| ( [2009] ) 
have a slope of around 0.8 (Reich, private communication). This 
is a strong indication that these fluctuations are mostly due to the 
magnetic field in M51. 

In Fig. [9] the RM fluctuations with an amplitude exceeding 
45radm~^, throughout the central r < 90", and 55radm~^ 
along the outer spiral arms cannot be explained by the noise. There 
are therefore around ten patches where RM changes sign over a 
distance of 1-2 kpc due to intrinsic fluctuations of magnetic field. 

The dispersion in RM is 15-20 rad m~^ measured in several 
regions in the inner spiral arms and 25-30 rad m~^ in the outer 
arms. After correction for the dispersion due to noise, the intrinsic 
dispersion is 1 1 =b 3 rad m~^ in the inner arms and 19 ± 5 rad m~^ 
in outer arms, hence constant within the errors. The distribution of 
rotation measures is shown in Fig.[To] along with the best-fit Gaus- 
sian, which has a mean RM of 10 db 1 rad m~^ and a dispersion of 
28 =b Iradm"^ 

Intrinsic fluctuations dominate the RM maps; even smoothed 
to linear scales of order 1 kpc, no large-scale pattern in rotation 
measure is apparent, in contrast to the clear large-scale spiral struc- 
ture in polarization angles. This result is quite surprising, as we 
would expect to see the components of the same field in the sky 
plane and along the line of sight in polarization angle and in Fara- 
day rotation, respectively. As polarization angles are not sensitive 
to field reversals, the observation of ordered pattern in angles does 
not demonstrate the existence of a regular (coherent) field. The spi- 
ral field seen in polarization angle could be anisotropic with many 
small-scale reversals, e.g. produced by strong shearing gas motions 
and compression, and hence would not contribute to Faraday rota- 
tion. Alternatively, the field may have significant components per- 
pendicular to the galaxy plane (due to loops, outflows etc.) which 
are mostly visible in Faraday rotation and hide the large-scale pat- 
tern. Such an underlying large-scale pattern indeed exists, as we 
discuss in the next section. 

The close alignment of the observed field lines along the 
CO arms and the lack of enhanced Faraday rotation in the polar- 
ized ridges can be understood if the turbulent magnetic field is 
anisotropic. An anisotropic turbulent field can produce strong po- 
larized emission, but not Faraday rotation. This picture is similar to 
that obtained for the effect of large-scale shocks on magnetic fields 
in the barred galaxies NGC 1097 and NGC 1365 ( |Beck et al.|2005l ) 
and will be investigated in detail in Sect. [7] 



10 A. Fletcher et al. 



5.1.1 The size of turbulent cells 

Here we derive a new method for estimating the size of turbulent 
cells in the ISM of external galaxies. 

The RM dispersion crRM,D observed within a beam of a linear 
diameter D is related to ctrm (Eq.|5]) as 



-1/2 



CTRM — CTRM 



D 



(2) 



where N — (D/d)^ is the number of turbulent cells within the 
beam area, assumed to be large. We confirmed the approximate 
scaling of (JiiM,D with using RM maps smoothed from 8^^ 
to 12^^ where the noise fluctuations are not dominant. Combina- 
tion with Eq. ^ allows us to estimate the least known quantity 
involved, the diameter of a turbulent cell (or twice the correlation 
scale of the turbulence): 



DaRM,D 



0.81(ne)Br(L)i/2 



2/3 



50 pc 



D 



600 pc 

Br 



20 /iG 



2/3 



-2/3 



/ cr-RM,D \ 
\ 15 rad m"^ / 



2/3 



{rie 



0.1 cm 



(3) 

-2/3 



1 kpc 



-1/3 



5.2 Faraday depolarization 

Faraday depolarization DP gives important information about the 
density of ionized gas, the strength of the regular and turbulent 
field components, and the typical length scale (or integral scale) 
of turbulent magnetic fields. DP is usually defined as the ra- 
tio of the degrees of polarization of the synchrotron emission at 
two wavelengths. This requires subtraction of the thermal emis- 
sion which is subject to major uncertainties (see Sect. |4]). Instead 
DP was computed, from the polarized intensities P, as DP = 



(Ph/Ph) X (zys/z^i)"^^", where ^s; 



- 1 . 1 is the synchrotron 



spectral index, assumed to be constant across the galaxy. Variations 
in asyn affect DP less severely than errors in the estimate of the 
thermal fraction of the total radio emission. The DP maps derived 
for A6 cm and A3 cm and between A20 cm and A6 cm are shown 
inFig.[TT] 

DP(6cm/3cm) (Fig. [TT^) is around unity (i.e. no Faraday 
depolarization) in most of the galaxy. Small patches with no- 
ticeable Faraday depolarization, where DP = 0.6-0.7 are gen- 
erally found in the spiral arms. There is no systematic connec- 
tion between the depolarization and the intensity of Ha emis- 
sion, indicating that variations in the thermal electron density 
are not the main source of the depolarization. The average value 
of DP(20cm/6cm) (Fig. ^) is 0.28, smaller by a factor of 
about 3 than DP (6 cm/3 cm). In the inner arms, DP(20cm/6cm) 
is lower than 0.2. Only in the outer regions of the disc does 
DP(20 cm/6 cm) increase to 0.5 or higher. 

Differential Faraday rotation within the emitting layer leads to 
depolarization which varies as a sin(x)/x, with x = 2 |RM^| A^, 
where RM^ is the intrinsic Faraday rotation measure within the 
emitting layer ( |Burn| 1 966| [Sokoloff e t al. 1 1 99 8| ) . With typical val- 
ues |RM^| = 50radm~ (Fig.|9|, we expect little depolarization 
(DP ^ 0.98) at A6.2cm and even less at A3. 6 cm. At A20.5cm, 
significant DP is expected for |RMi | > 30 rad m~^. Furthermore, 
lines of zero polarization ('canals') are expected along level lines 
with |RMi| = 37.5nradm~^ (with integer n 7^ 0) (Shukurov 



|& Berkhuijsenl[2003l [Fletcher & Shukurov|[2008] ), but not a sin- 
gle 'canal' is observed in the A20.5cm polarized intensity map. 



This suggests that the average \RMi\ at A20.5cm is significantly 
smaller than 37.5 rad m~^. Since the depolarization is relatively 
strong, it must be due to a different mechanism, e.g., Faraday dis- 
persion. 

Internal Faraday dispersion by turbulence in the magneto- 
ionic interstellar medium is the probable source of strong depolar- 
ization at longer wavelengths, producing the degree of polarization 
given by ( |Burn|1966| [Sokoloff et al.|1998l ): 



■Po- 



1 -exp(-25} 
25' 



(4) 



where S — cr^^A^ and the maximum degree of polarization is 
po ~ 0.7. Here ctrm is the dispersion of intrinsic rotation measure 
RMi within the volume traced by the telescope beam. 



fTRM =0.8l(ne)Br(L(i) 



1/2 



(5) 



where (rie) is the average thermal electron density along the line of 
sight (in cm~^), Br the strength of the component of the random 
field along the line of sight (in /iG), L the total path-length through 
the ionized gas (in pc), d the size (diameter) of a turbulent cell (in 
pc). Reasonable values for the thermal disc of (rie) = 0.1 cm~^, 
L = 800 pc (estimated for 2.4 < r < 4.8 kpc by Berkhu ijsen| 
|eraLiri997 ), Br = 20 /xG (Sect. |4}, d = 50 pc (see Sect. ^A) 
yield ctrm ~ 300 rad m~^ resulting in p/po ~ 0.002, a very 
strong depolarizing effect. The resulting Faraday depolarization is 
DP = 0.002 between A20.5 cm and A6.2 cm, much smaller than 
the observed depolarization of DP = 0.2. 

We conclude that the A20 cm polarized emission must origi- 
nate from a layer at a greater height in the galaxy than the bulk of 
the A3 cm and A6 cm polarized emission, as the A20 cm disc con- 
tribution to the polarized signal that we observe in Fig. p] must be 
severely depolarized. Two possibilities, that are not mutually exclu- 
sive, are that (i) the scale height of the A20 cm synchrotron disc is 
greater than the scale height of the thermal electron disc, or (ii) the 
A20 cm polarized emission is produced in a synchrotron halo. In ei- 
ther case, since Faraday rotation is observed at wavelengths around 
A20cm ( [Horellou et al.|1992[|Heald et al.|2009) , with about 1/5 
the amplitude as between A3 cm and A6 cm (see below), thermal 
electrons and a magnetic field are required in the halo0 

In contrast, at A3. 6 cm ctrm ~ 300 rad m~^ produces virtu- 
ally no depolarization with p/po ~ 0.85, so we expect that the 
galaxy is transparent (or Faraday thin) at A3. 6 cm. At A6.2cm 
we have p/po ~ 0.3, moderate depolarization, with DP = 0.2 
between A6.2cm and A3. 6 cm. This rough estimate for DP is 
the same order of magnitude as the DP shown in Fig. [TTJa), al- 
beit about three times lower than the typical observed value of 
DP ^ 0.9, indicating that our value of ctrm is a slight overes- 
timate. (The uncertainty in the adopted values of (rie), L and Br 
can easily explain the discrepancy: for ctrm ~ 200 rad m~^ we 
have DP = 0.6, much closer to the observed values.) Since Fara- 
day effects are so small at A3. 6 cm. Fig. [TTJa) strongly indicates 
that the disc is also Faraday thin at A6.2 cm. 



^ The terminology used here is unavoidably imprecise. Instead of a disc 
and halo we could just as easily refer to a thin and thick disc. We are unable 
to say anything about the geometry (e.g. flat or spheroidal) of the two layers 
using our data, only that there must be two layers producing the Faraday ro- 
tation observed in M51. In Section[6]we look at the magnetic field structure 
in the two layers in more detail. 
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Figure 11. Depolarization between (a) A6cm and A3 cm and (b) between A20cm and A6cm, both at 1^" resolution. Also shown are contours of Ha 
emission jGreenawalt et al.|1998 ) at the same resolution, plotted at 4, 8, 16, 32, 74% of the map maximum. 



6 REGULAR MAGNETIC FIELD STRUCTURE 
6.1 The method 

The map of Faraday rotation discussed in Sect. |5.1| clearly shows 
the effects of strong magnetic field fluctuations on scales of 300 pc 
to ^ 1 kpc; from this map one could expect the magnetic field to be 
rather chaotic and disordered. However, the observed polarization 
angles, shown in Figs. [T] and [2] suggest an underlying spiral pat- 
tern to the magnetic field on scales > 1 kpc even when corrected 
for Faraday rotation (Fig.|4]). If the large-scale spiral pattern in the 
magnetic field is due to a regular field component, such as might be 
expected due to mean-field dynamo action, we can expect to find a 
signature of such a field in the Faraday rotation signal at relevant 
scales. If no such signal can be uncovered then we may conclude 
that the data do not support the presence of a mean field in M5 1 and 
that the spiral patterns seen in Figs.[T]and[2]are being imprinted on a 
purely random magnetic field through e.g. large-scale compression 
in the gas flow. 

\ shows a Faraday rotation map using our A3 cm and 
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A6cm data smoothed to 30^^ resolution: large-scale structure in 
the rotation measure distribution across the disc is now visible (cf. 
Fig. The rotation measures between A18cm and A20cm also 
show large-scale structure ( [Horellou et al.|1992 Fig. 9), but with 
a markedly different pattern. The complex regular magnetic field 
structure producing the observed rotation measure patterns cannot 
be reliably determined from an intuitive analysis of these maps. For 
example, a simple axisymmetric or bisymmetric field would pro- 
duce a single or double period variation with azimuth, whereas the 
actual pattern is clearly more complicated. A particular difficulty 



is that the pattern in the observed Faraday rotation is very different 
between the pairs of short and long wavelengths, indicating that 
two Faraday-active layers may be present. 

In order to look for the signature of a regular magnetic field 
in the multi-frequency polarization maps we applied a Fourier fil- 
ter to the 15^^ resolution Stokes Q and U maps at A A3, 6, 20 cm 
to remove the signal on scales < 30^^ ^ 1.1 kpc. We also used 
the A18 cm data of [Horellou et al.| ( [l992| ) which has a resolution of 

resolution and was not filtered. Then maps of polarization an- 
gle were constructed at each wavelength. These were subsequently 
averaged in sectors with an opening angle of 20° and radial ranges 
of 2.4-3.6 kpc, 3.6-4.8 kpc, 4.8-6.0 kpc and 6.0-7.2 kpc (see 
Fig.[T4|. We estimated that the minimum systematic errors in po- 
larization angle arising from this method are about 4° at A A3, 6 cm 
and 10° at AA18, 20 cm, with the main source of error being un- 
certainties arising from Faraday rotation by the random magnetic 
field (see Ruzmaikin et al. 1990 Sect. 2). We set these as mini- 
mum errors in the average polarization angles, otherwise using the 
standard deviation in the sector. 

We applied a method that seeks to find statistically good fits to 
the polarization angles using a superposition of azimuthal magnetic 
field modes exp(im(/)) with integer m, where (j) is the azimuthal 
angle in the galaxy's plane measured anti-clockwise from the north 
end of the major axis. A three-dimensional model of the regular 
magnetic field is fitted to the observations of polarization angles 
at several wavelengths simultaneously. The polarization angle af- 
fected by Faraday rotation is given byV^ = V^o + RMA^ + RMfgA^ 
where the intrinsic angle of polarized emission is ^o, RM is the 
Faraday rotation caused by the magneto-ionic medium of M5 1 and 
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Figure 12. Rotation measures between A A3, 6 cm, at 30^^ resolution, over- 
laid with contours of mid-infrared 15 /im emission ( Sauvage et al.|199"6) 
plotted at 1, 2, 4, 8, 16, 32% of the maximum value. Data were only used 
where the signal-to-noise ratio in polarized intensity exceeds five. 



RMfg is foreground Faraday rotation arising in the Milky Way. 
The method of modelling is described in detail in |Berkhuijsen et| 
|ar| ( |1997| ) and |Fletcher et al.| ( |2004|) and has been succe ssfully ap- 
plied to normal spiral galaxies ( [Berkhuij sen et al.| 1 997 [ [Fletcher et 



arf2004 , Tabatabaei et al. 2008| and barred galaxies ( Moss et al. 



2001||Beck et al.|2 005 ). We have derived a new model of the regu- 



lar magnetic field in M5 1 as our combined VLA -i- Effelsberg maps 
at A A3, 6 cm are a significant improvement in both sensitivity and 
resolution over the maps used by [Berkhuij sen et aL] ( |1997| ). 

The fitted parameters of the regular magnetic fields are given 
in Appendix |A| These fit parameters can be used to reconstruct 
the global magnetic structure in M51. In order to obtain statisti- 
cally good fits to the observed data we needed plane-parallel field 
components only (no vertical fields). No satisfactory fits could be 
found, using various combinations of the azimuthal modes m=0, 1, 
2 for the horizontal field and m=0,l for the vertical field, for the 
observations at all 4 wavelengths for a single layer. Therefore at 
least two separate regions of Faraday rotation are required. This 
is because the patterns of polarization angle and Faraday rotation 
at A A3, 6 cm and A A 18, 20 cm are very different: at A 18 cm and 
A20 cm the disc emission is heavily depolarized by Faraday dis- 
persion (see Sect. [5.2| ), so we only see polarized emission at these 
wavelengths from the top of the disc. A similar requirement for two 
Faraday rotating layers in M5 1 was also found by [Berkhuij sen et[ 
|ar||1997). 

To describe the two layers in the model, the Faraday rotation 
from M51 is split into two components, arising from a disc and 
halo, RM = ^(^)RM(^)+^(^)RM(^\ where ^(^^ and^^^^ are 
parameters that allow us to model how much of the disc and halo 




Figure 13. Geometry of the disc and halo layers. The thermal disc and 
halo have thickness /ith and Z, and electron densities n^^d and n^^h- 
The scale-height of the region emitting polarized synchrotron radiation at 
AA18, 20 cm is /i2o and the observed A20cm polarized emission comes 
from the layer /ith < z < h2o, due to Faraday depolarization. The regular 
magnetic field has two layers: the disc B^, extending to the same height as 
h2o, which is both an emitting and Faraday rotating layer, and the halo B]^ 
where only Faraday rotation occurs. 



are visible in polarized emission at a given wavelength. We use this 
decomposition of the RM into disc and halo contributions to take 
into account the depolarizatio n of the A20 cm emission from the 

' ]by setting 



thermal disc discussed in Sect. |5.2| by setting ^''^^ = and = 
1 at AA18, 20 cm and ^^^^ = 1 and ^^^^ = 1 at AA3, 6 cm. In 
other words, at AA18,20cm the polarized emission is produced 
in a thin layer that lies above the thermal disc (see Sect. \5.2\ and 
has the same regular magnetic field configuration that produces the 
A A3, 6 cm polarized emission. The A A3, 6 cm emission is Faraday 
rotated in the thermal disc and the halo whereas the A A 18, 20 cm 
emission is only Faraday rotated in the halo: see Fig. [13] 



6.2 Results 

The resulting regular magnetic field structure is shown in Fig. [14] 
The regular field in the disc is best described by a superposition of 
m=0 and m=2 horizontal azimuthal modes and has a radial com- 
ponent directed outwards from the galaxy centre, whereas the halo 
field has a strong m = 1 horizontal azimuthal mode and is directed 
inwards in the north, opposite to the direction of the disc field, and 
outwards in the south, same as the disc field. Oppositely directed 
components of the field in the disc and halo were also found by 
jBerkhuijsen et al.| ( |1997| ), however our new observations place the 
strong m = 1 mode, and the resulting strong asymmetry in the 
magnetic field, in the halo rather than the disc. We believe that the 
difference in fitted fields arises from the higher quality of the new 
A A3, 6 cm data, as discussed above. Heald et al.| ( [2009) derived ro- 
tation measure maps from multi-channel WSRT data at A22 cm us- 
ing the RM- synthesis technique that qualitatively show the pattern 
expected from a m = 1 magnetic field. Since depolarization due to 
Faraday dispersion cannot be removed by RM-synthesis, this gives 
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Figure 14. The regular magnetic field derived from fitting a model to the observed polarization angles at A A3, 6, 18, 20 cm, with the length of the magnetic 
field vectors proportional to the field strength, overlaid on the same optical image as in Fig. (a) The regular magnetic field in the galactic disc, (b) The 
regular magnetic field in the galactic halo. Ring boundaries are at 2.4, 3.6, 4.8, 6.0 kpc and 7.2 kpc and all sectors have an opening angle of 20°. The major 
axis is indicated: the midpoints of these two sectors correspond to (/) = 0° and (f) = 180° respectively. 



a second, independent, indication that the halo of M51 hosts an 
m = 1 regular magnetic field. 

In the ring 2.4 < r < 3.6kpca weak m = mode is required 
to fit the data. This arises due to a sharp change in the observed po- 
larization angles at AA18, 20 cm between = 120° and = 140°. 
Even with two halo modes we cannot capture the rapid change in 
angles and had to exclude the A20 cm data at = 140°: the alter- 
native would be to add an extra mode, with three new parameters, 
to model one data point. 

The process by which two different regular magnetic field pat- 
terns in two layers of the same galaxy are produced is not clear and 
is beyond the scope of this paper. We only offer some speculative 
suggestions: (i) both fields might be generated by mean field dy- 
namos but operating in different regimes, with the interaction of 
M51 with NGC 5195 driving a m = 1 mode in the halo; (ii) the 
halo field could be a relic of the magnetic field present in the ten- 
uous intergalactic medium from which the galaxy formed; (iii) as 
the disc field is advected into the halo it can be modified by the 
halo velocity field into the m = 1 pattern. All of these possibilities 
will require careful modelling to determine their applicability to the 
problem. 

In the disc of M5 1 , in the four rings used in our model the m — 
azimuthal field component is 1-2 times the strength of the m = 2 
mode (Table [aT| . While the strength of the m = 2 mode remains 
approximately constant between the rings, the m = mode is of 
similar strength in the inner ring, but is much stronger in the other 
three rings. 



The r.m.s. regular field strength B in each ring can be deter- 
mined by integrating the fitted modes over azimuth (Table [aTJ via 

88radm-2 V0.11cm-3y Vlkpcy ^ ' 
where 

with Rr and i?^ given by Eq. ( [Al| ). |Berkhuijsen et al.| { [1997] ) es- 
timated that (ne,d) = 0.11 cm~ and hth — 400 pc in the ra- 
dial rang^ 2.4 < r < 4.8 kpc and for 4.8 < r < 7.2 kpc 
they estimated (ne,d) = 0.06 cm~^ and hth — 600 pc. The r.m.s. 
strengths of the large-scale magnetic field, using these parameters, 
are5 = 1.4±0.1/iG, B = 1.7±0.5/iG, B = 2.7±1.0/iG 
and 5 = 2.8 ± 0.1 /xG in the rings 2.4 < r < 3.6 kpc, 
3.6 < r < 4.8 kpc, 4.8 < r < 6.0 kpc and 6.0 < r < 7.2 kpc 
respectively. These are a factor of 4 lower than the strength of the 
ordered field derived from the equipartition assumption (Sect. |4.2| . 
The equipartition field strength is based on the observed polarized 
intensity and an anisotropic random magnetic field can contribute 
to the polarized signal (see Sect. |4.2| ) but will not produce any sys- 
tematic pattern in polarization angles at different frequencies. This 
is probably the reason for the discrepancy: most of the polarized 

^ [ Berkhuijsen et al.Hl997t adopted a distance to M51 of 8 Mpc: we have 
rescaled their radii to our distance of 7.6 Mpc. 
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radio emission in M51 does not trace a mean magnetic field, only 
the modelled large-scale pattern in Faraday rotation does. 

The average pitch angle of the m = mode is —20° with 
little variation in radius between the rings. This means that the spi- 
ral structure of the regular field is coherent over the whole galaxy. 
The weaker m = 2 modes produce an azimuthal variation in pitch 
angle of about 15° in the inner ring and 5° in the other rings. This 
variation of the pitch angle of the regular magnetic field with az- 
imuth is much lower than the variation of the pitch angle of the 
ordered magnetic field observed in polarization: [Patrikeev et al.| 
P006| l showed that the orientation of Faraday rotation corrected 
polarization angles change by about 30°. Thus the anisotropic ran- 
dom magnetic field, that we believe produces most of the polarized 
emission in M51, has a stronger azimuthal variation in its orienta- 
tion: this is to be expected if the anisotropy arises from some pe- 
riodic mechanism such as compression in spiral arms or localised 
enhanced shear. 

The bisymmetric halo field has a much larger pitch angle of 
about —50° in the inner three rings. If this field is generated by a 
mean field dynamo the high pitch angles may be an indication that 
differential rotation in the halo is weak and an -dynamo action 
is significant. (In an -dynamo, shear due to differential rotation 
is negligible and so |Br| ~ l^^l) The average thermal electron 
density and size of the halo in this galaxy are unknown and we 
have no specific constraints to apply. Taking reference values of 
(ne,h) 0.01 cm"^ Z = 5kpc and (ne,h) 0.06 cm"^ Z 
3.3 kpc for the radial ranges 2.4 ^ r ^ 4.8 kpc and 4.8 ^ r ^ 
7.2 kpc respectively, where the densities are one tenth of the disc 
density as in the Milky Way and Z is as estimated by |Berkhuijsen et| 
[ar] ( |1997| , the r.m.s. strength of the halo field (note that the m = 1 
mode means there are two values of the azimuth in each ring where 
the field is zero) is 



B : 



Rh 



{ne,h) 
^radm~2 V0-01cm~ 



1 kpc 



where Rh is the average amplitude of the halo mode. The fitted 
amplitudes of the halo field given in Table |Al| give r.m.s. regular 
field strengths in the halo of 1.3 ± 0.3, 1.2 ± 0.4, 2.2 ± 1.0 and 
1.6 it 0.6 /iG in the four rings, with increasing radius. 

We are confident that the azimuthal modes and pitch angles, 
fitted to our data, for the regular magnetic field in the disc and halo 
are robust as we have carried out extensive checks and searches 
of the parameter space. For example, if the AA18, 20 cm data is 
ignored and a disc only model used the A A3, 6 cm data produce 
a very similar fitted field to that given in Table |A1| However, we 
cannot consider the mode amplitudes to be reliable other than to 
reflect the relative strengths of the m = and m — 2 modes in the 
disc. 

A field reversal between the regular fields in the disc and the 
inner halo has also been suggested in the Milky Way ('Sun et al.| 
|2008| ), perhaps similar to the northern half of M51 (Fig.[T4jr 

It is probable that our model for the vertical structure of M5 1 
(Fig. [13]) is too simple and that this has lead to too much Faraday 
rotation being put into the halo field at the expense of the disc. In 
particular we do not allow for Faraday rotation from the thin emit- 
ting part of the A20 cm disc nor for emission at any wavelength 
from the halo. Furthermore, each azimuthal mode is assumed to 
have an azimuth-independent intrinsic pitch angle. Adding more 
parameters to this model, given the limited number of sectors and 
frequencies that we can use, will not likely resolve this question. A 
more productive approach will be to develop a new model that also 
takes into account depolarizing effects directly and whose outputs 



are statistically compared directly to the individual maps, including 
the unpolarized emission (perhaps using the maps of Stokes param- 
eters themselves). 



7 ARM-INTERARM CONTRASTS 

We wish to investigate the effect of magnetic field compression in 
the large-scale shocks associated with the spiral arms. In partic- 
ular, how the regular and random magnetic field components are 
changed by the shocks and whether shock compression of isotropic 
random magnetic fields can produce enough anisotropic field to ac- 
count for most of the polarized emission as inferred in Sections [5^ 
and [6| We have carefully examined the azimuthal variations in 
the various maps at different radii: these emission profiles clearly 
demonstrate that there is not a simple azimuthal behaviour of any 
of the measured quantities. One cannot identify "typical" arm to 
interarm contrasts. So we have used a mask to separate arm (pixels 
in the mask set to 1) and interarm (mask pixels set to —1) regions 
in each of the maps and hence calculate the average contrast over a 
wide radial range. 

We combined the CO map of [Heifer et aT] ( [20031 ) with the 
HI map of Rots et al.| ( |1990| ) to produce a map of the total neutral 
gas density at 8^^ resolution, assuming a constant conversion fac- 
tor A^H2 = 1.9 X 10^*^/00 cm~^. The arm-interarm mask was 
determined by making a wavelet transform, using the Mexican hat 
wavelet with a linear scale of about 1 kpc, of this map. This scale 
was selected by examining a range of transform maps: a 1 kpc scale 
wavelet produces transform coefficients that are continuously pos- 
itive along spiral arms and continuously negative in interarm re- 
gions. Also, 1 kpc seems reasonable as a typical width of the spiral 
arms in M5 1 . This mask was used to separate the arm and interarm 
components of each of the maps listed in Table |2] over the radial 
range 1.6 kpc < r < 4.8 kpc, and then the average arm and inter- 
arm values were determined. 

In addition to the data in Table[2|we also separated the rotation 
measures shown in Fig. [9] into arm and interarm components. We 
calculated the average magnitude (|RM|) and standard deviation 
CTRM of RM and found that (|RM|) ^ 22radm~^ andcTRM ^ 39 
in the interarm and (|RM|) ^15 rad and ctrm ~ 29 in the 
arms. This may indicate that the regular magnetic field is stronger in 
the interarm regions than in the arms. However, the interpretation 
is difficult as the RM distribution depends on the signal-to-noise 
ratio, which tends to be higher in the arms. 

The contrast in the neutral gas column density (2H2 + HI) is 
compatible with what might be expected from compression by a 
strong adiabatic shock, 

where the superscripts and refer to downstream and up- 
stream of the shock front. We note that the scale height of the 
gas layer h is not expected to be much affected by the spiral pat- 
tern ( |Shukurov| l998); Hi observations in t he Milky Way sugge st 
^(d) j^w ^ i_i_5 in the outer Milky Way |Kalberla etal.|2007l. 



7.1 General considerations 

Explaining the arm-interarm contrast in the observed radio inten- 
sity is a long-standing problem. |Mathewson et al.| ( |1972| ) argued 
that their A20 cm observations of M51 are consistent with the den- 
sity wave theory of [Roberts & Yuan| ( |1970| ). They assumed that 
both cosmic ray number density and the tangential magnetic field 
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Table 2. Arm-interarm contrasts in observed quantities. Radio intensities at A6 cm are shown. The data has been smoothed to 1^" resolution and the arm and 
interarm regions identified using a mask derived from the wavelet transform of the combined CO and HI map (see the text for details). The inner arm data are 
for a resolution of A" at a position of strong contrast in polarized radio emission. 



Quantity 


Units 


Average 1.6 kpc < r < 4.8 kpc 


Inner arms 0.8 


kpc < r < 1.6 kpc 






Arm 


Interarm Arm/interarm ratio 


Arm 


Interarm 


Arm/interarm ratio 


Neutral gas column density 


102iHcm-2 


18.0 


4.0 4.5 


200 


40 


5 


Total radio intensity 


mJy/beam 


1.1 


0.5 2.2 


0.6 


0.12 


5 


Polarized radio intensity 


mJy/beam 


0.1 


0.1 1.0 


0.07 


< 0.01 


^ 7 



increase in proportion to the gas density at the spiral shock. The 
resulting arm-interam contrast in radio intensity, after taking into 
account of the telescope beam width, is expected to be of order 10 
or more. 

|Tilanus et al.| (|1988), using observations at a higher resolution, 
found that the shape of cross-sectional profiles across the radio in- 
tensity arms is not compatible with the density wave theory and 
concluded that the synchrotron emitting interstellar medium is not 
compressed by shocks and decouples from the molecular clouds as 
it traverses the arms. Thus there is clearly a discrepancy between 
the physically appealing theory of [Roberts & Yuan| ( |1970| ) and ob- 
serv ations (see also [Condon] 1992| p. 590). 

jMouschovias et al. ( 191 A) and jMouschovias et aL] ( [2009| ) sug- 
gest that only a moderate increase in synchrotron emission in spi- 
ral arms is expected due to the Parker instability: rather than being 
strongly compressed the regular magnetic field rises out of the disc, 
in loops with a scale of 500-1000 pc. However, the substantial ran- 
dom component of the magnetic field in M5 1 may supress the insta- 
bility or reduce it to a simple uniform buoyancy ( Kim & Ryu 2001). 
Furthermore we do not observe the periodic pattern of enhanced 
Faraday rotation along the spiral arms that would be expected from 
the vertical magnetic fields at the loop footpoints (Fig. [9]). 

In this section we reconsider this question with additional em- 
phasis on the polarized intensity. In Sect. [7.2[ we consider how 
shock compression affects the emission of cosmic-ray electrons at 
di fixed frequency. The effect of the modest observational resolution 
on the arm-interarm contrasts is discussed in Sect. 17.31 In Sect. 17.41 
we consider the case of compression of an isotropic random mag- 
netic field and conclude that this may be the dominant origin of the 
arm-interam contrast in radio intensity only in the inner galaxy. Fi- 
nally, in Sect. [73] we show how the decompression of an isotropic 
random magnetic field as it leaves the spiral arm affects the arm- 
interarm contrast. 



7.2 Cosmic rays in compressed gas 

Assuming that magnetic field is parallel to the shock front of the 
spiral density wave and is frozen into the gas, its strength increases 
in proportion to the gas density, B oc p, as appropriate for one- 
dimensional compression. The ultrarelativistic gas of cosmic rays, 
whose speed of sound is c/v^, with c the speed of light, is not 
compressed in the arms. However, compression of magnetic field 
will affect the cosmic rays (including the electr on component) be- 
cause p']_/B ^ const is an adiabatic invariant ( [Rybicki & Light- 
man ][T979'), where p± — jiTieC is the component of the particle 
momentum perpendicular to the magnetic field and 7 is the particle 
Lorentz factor. More precisely, only the part of the Lorentz factor 
related to the particle velocity perpendicular to the magnetic field 
should be included, but we ignore this detail for a rough estimate. 
In terms of the Larmor radius r^, one can write p± = eBrs/c, 



with e the electron charge, to obtain another form of the adiabatic 
invariant, Br% = const, i.e. magnetic flux through the electron's 
orbit remains constant. For B oc p, we then obtain rs oc p~^^^ 
and p± oc p^^^, or 7 oc p^^^. 

Thus, compression of magnetic field leads to an increase in 
the Lorentz factor of the cosmic-ray electrons, 7 oc p^^^. If the 
initial range of the Lorentz factors is 7^-^^ ^ 7 ^ 7mix, compres- 



sion transforms it into 7^/^ ^ 7 ^ 7max such that 7min/7i 
(^(d)/^(u))l/2^ Of course, the total number density of cosmic-ray 
particles does not change because the cosmic-ray gas is not com- 
pressed. Adopting a power-law spectrum of cosmic-ray electrons, 

d7 = K^^~^ d7 , 



where d^y is the number of relativistic electrons per unit volume 
in the range (7, 7 + d7), the total number density of the particles 
follows as J^^^^ d7 ^ 7m7n ^7 = const, where we have 

assumed that the energy spectrum is broad enough to have 7max ^ 
7min and s > 1. However, the energy of each cosmic-ray particle 
increases as p^^^: the energy shifts along the energy (or 7) axis, 
and 

i.e., the number density of particles with given 7 increases with p. 

Now we can estimate the effect of compression on the syn- 
chrotron emissivity observed at a fixed frequency, v — const and 
fixed frequency interval — const. Denoting ei(y) the arm- 
interarm contrast in liv) we have 



OC p 



since / oc K^B^'^^^^^ 



^ 7n 



-i^(s+i)/2^ Note that the Lorentz 



factor of the electrons which radiate at a fixed frequency ly 
j^B = const reduces as B increases, 7(z^) oc this also 

leads to an increase in the number of cosmic-ray electrons radiat- 
ing at a given frequency after compression since rij oc 7"* and 
s > 0. 

With s ^ 3 and the arm-interarm density contrast of about 
four, the number of cosmic-ray electrons with a given 7 is propor- 
tional to Kj oc p*^*"^-*/^ oc p, and the synchrotron intensity in 
the spiral arms would then be 50-100 times stronger than between 
the arms. There are other reasons to expect enhanced synchrotron 
emission in the arms as supernova remnants, sites of cosmic ray 
acceleration, are localized in the arms. However, Table |2] clearly 
shows that such an enormous contrast is not observed. 

It is plausible that cosmic rays are rather uniformly distributed 
in galactic discs and only weakly perturbed by the spiral arms (Sec- 
tion 3.10 of jBerezinskii et al.|1990| ). During their lifetime within 
the galaxy, r ^ 3 x 10 yr, the cosmic ray particles become 
well mixed over distances of order (Dr)^/^ ^ 2-3 kpc, where 
D 4 X 10^^ cm^ s~^ is the cosmic ray diffusivity. This scale ex- 
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ceeds the width of spiral arms, so diffusion can significantly reduce 
the arm-interarm contrast, but it can hardly result in the almost uni- 
form distribution of cosmic rays which would explain the low con- 
trasts in Table|2] Anyway, even assuming that the cosmic ray inten- 
sity is the same within the ams and between them, the compression 
of magnetic field by a factor of four would result in an enhance- 
ment of the synchrotron emissivity by a factor 16 for s = 3, and 
this is already larger than what is observed. Note that although the 
contrast in total radio emission given in Table |2] includes the ther- 
mal radio emission, this is concentrated in the arms as can be seen 
in Ha images fGre enawalt et al.|1998| ) and so the actual contrast in 
synchrotron emissivity will be lower than the contrast in total radio 
emissivity. 

7.3 Telescope resolution and the width of the compressed 
region 

One factor that can help to explain the lower than expected arm- 
interam contrast in total radio emission (which for simplicity we 
assume to be all synchrotron emission) is the possibility that the 
compressed region is narrow compared to the width of the spiral 
arms. Given the complex processes that take place as gas passes 
through the arms, such as the formation of molecular clouds and 
star formation, it is plausible that shock compression is followed by 
a decompression before the gas leaves the arms. Then, if the com- 
pressed region is narrower than the telescope beam the observed 
contrast between the arm and interarm will be reduced. We can es- 
timate the width of the compressed region that is compatible with 
the observations in Table [2]using 



where ej 5 is the observed arm-interam contrast in total radio 
emission in the inner spiral arms (which we assume to be all syn- 
chrotron emission to make a conservative estimate), :^ 50 is the 
expected contrast due to compressive effects in an adiabatic shock, 
D ^ 150 pc is the beam width at our highest resolution and w is the 
width of the compressed region. Using the values of the parameters 
just quoted we obtain w ^ 15 pc. 

|Patrikeev et al^ ( [2006 )) showed that the ridge of strongest polar- 
ized emission (tracing the peak of the compressed magnetic field) 
is generally shifted upstream of the ridge of strongest CO emission 
(tracing the highest neutral gas density) in M51. We can assume 
that this shift is due to a spiral- shock triggering the formation of 
molecular clouds. The dense clouds will fill only a small fraction 
of the volume occupied by the ISM and may become decoupled 
from the magnetic field threading the diffuse ISM during their for- 
mation, as originally suggested in the case of M51 by |Tilanus et| 
laTI (19881 ), by [Beck et al.| ((2005) for barred galaxies and with a 
plausible mechanism for the separation outlined in F letcher et al.| 
( |2009| ). We expect that the expansion of the compressed magnetic 
field will begin once the clouds have formed. Thus after a time 
Tc ~ 10^ yr the ridge in strong radio emission will begin to de- 
cay. If we estimate the magnetic field strength to be B 20 /xG 
(Sect. [4| then an Alfven wave will propagate over the compressed 
distance w = 15pc if the density, of the compressed diffuse gas 
in which the clouds are embedded, is n 8 cm~^. This density is 
not implausible if the upstream diffuse gas density is n :^ 2 cm~^. 

So one possible explanation for the observed low arm-interam 
contrast in total radio emission, compared to that expected from a 
simple consideration of cosmic ray energies, is that the compressed 
region is narrower than the beam. In this case we estimate that the 



ridge in compressed magnetic field, along the upstream edge of the 
spiral arm, will be a few tens of parsecs wide. This possibility can 
be tested using higher resolution observations. The cosmic rays are 
not the only source of the arm-interam contrast in radio emission 
though: next we consider the effect of a large-scale shock on the 
regular and random components of the magnetic field. 

7.4 Compression of a partially ordered magnetic field 

We consider compression of a partially ordered magnetic field by 
the spiral shock. We assume that the random magnetic field up- 
stream of the shock is statistically isotropic but one-dimensional 
compression makes it anisotropic and so it contributes to the po- 
larized radio emission but not to Faraday rotation. We assume that 
both the spiral arms and the field lines of the large-scale magnetic 
field are logarithmic spirals with the pitch angles pa and pb, respec- 
tively. An acceptable estimate is Pa ~ P6 ~ —20°, where negative 
values correspond to a trailing spiral. 

The formulae required to calculate the compression of a mag- 
netic field with both regular and isotropic random components and 
the associated total and polarized synchrotron emission and Fara- 
day rotation are derived in Appendix [B] Using these equations we 
now consider two cases that encompass the range of observed con- 
trasts in polarized emission in M51: the inner spiral arms, where 
there is a strong arm-interarm contrast in polarization of at least 
epi ^4, and the average contrast of epj ^ 1 which is more rep- 
resentative of the greater radii. The upstream ratio of random to 
regular magnetic field strength b/B and the increase in gas density 
Cn = 4, in regular field cb and in random field Cb fields due to the 
shock are fixed and then the consequential increases in nonthermal 
6/ and polarized epi emission across the shock are calculated. 

Firstly the inner arms. Here we can obtain a reasonable match 
with the maximum observed arm-interarm contrasts ej ^ 5 (as- 
sumed to be mostly synchrotron emission) and epi 7 (a lower 
limit due to weak interarm intensity at 4'' resolution). With the 
parameters = 1, i.e. the regular field is not increased by the 
shock — this can be justified, for example, if the regular magnetic 
field becomes detached from dense clouds as they form ( [Beck et aL] 
|2005| ) — and Cb = 2.7 we obtain expected arm-interarm contrasts 
in synchrotron emission of 3.3 and polarized emission of 4.8. If 
> 1 then the expected increase in polarization becomes much 
larger. For example, for = = 2.7 we obtain epi = 11: 
since we only have a lower limit on the observed e pi we cannot 
rule out this possibility. We conclude that in the innermost spiral 
arms (r < 1.6 kpc) anisotropic random magnetic fields produced 
by compression of the interstellar gas in the spiral arms can account 
for the observed increases in total and polarized radio emission. 

Now for the average arm-interarm contrasts. Here, since the 
average contrast in polarization is epj ^ 1, this requires no in- 
crease in regular magnetic field in the arm and also no anisotropy of 
the random magnetic field, but simultaneously we require a small 
increase in synchrotron emission e/ = 2.2. This is not possible 
in our model; it is also unlikely to happen when galaxies contain 
strong spiral density waves. The closest we can come is to set 
= 1 coupled with a weak upstream random field / — 0.2 
and strong compression of the random field tb = 4.5 to produce 
the required contrast in total emission. Then we can obtain arm- 
interarm contrasts of 2.3 in total synchrotron emission and 1.7 in 
polarization. However, }? jB'^ — 0.2 is strongly contradicted by 
the equipartition estimates of the various field strengths in Sect. |4] 
which likely overestimate the strength of the regular field as dis- 
cussed in Sect. l5.2l and[6l 
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To summarise this subsection: enhancement of the regular and 
random magnetic field components parallel to a large-scale spiral 
shock can partly account for the observed arm-interam contrasts in 
radio emission but no single set of parameters is compatible with 
the full range of the observations. 

7.5 Decompression of an isotropic magnetic field 

Finally, we consider the decompression of an isotropic field as the 
magnetized gas flows from a high density to a low density region. 
This decompression also leads to the generation of an anisotropic 
random magnetic field, as the field component parallel to lines 
of constant density increases while the perpendicular component 
is unaffected, similar to the compressive case. But decompres- 
sion will also help to lower the arm-interam contrast in nonther- 
mal emission, thus alleviating one of the difficulties encountered in 
the previous subsection, if the random field in the arms is (partly) 
isotropized: this can readily occur due to turbulence being driven by 
star formation activity and the expansion of supernova remnants. 

First consider the case of straightforward compression. Let 
us define the x-direction as perpendicular to the shock and the y- 
direction as parallel: so only ^-components of the magnetic field 
are affected by the shock. Thus 

3 

and 

j2{d). _ 2n2{u). 1 







2i2{u) 



where {u) and (d) refer to upstream and downstream and b is the 
total random field. Now the plane-of-sky component of the random 
field is h\ — hi + so we obtain 

for an adiabatic shock, where Cp = 4. So we expect adiabatic com- 
pression in a simple, plane parallel shock to produce a contrast in 
nonthermal emission e/ 8, in the case that the cosmic rays are 
smoothly distributed. This is a simplified version of part of the cal- 
culation described above and in detail in Appendix [B] 

Now consider the case of decompression. We follow a similar 
calculation but here 



U2(u) 



g2 \''y 



which leads to 

In this case the arm is the upstream region so the expected arm- 
interam contrast in nonthermal emission is ej 2 which is close 
to the average contrast observed (Table [2]). 

7.6 Summary: compression and arm-interarm contrasts 

In Sect. |5.2| and[6]we have shown that much of the polarized emis- 
sion detected across the disc of M51 must come from anisotropic 
random fields. Combined with the problem described in Sect. |7.4| 
above, of how isotropic random fields can be compressed in a spi- 
ral shock but not produce an increase in polarized emission, we are 
led to the view that anisotropic random fields are already present in 
interarm regions, perhaps as a result of enhanced localised shear or 



decompression. [Patrikeev et ar]p006| ) showed that the orientation 
of the magnetic field in M51 varies in azimuth by ±15° and in the 
interarm generally has a different pitch angle to the CO spiral arm, 
only becoming well aligned with the CO arm at its location. In this 
case compression of the already anisotropic field in the spiral shock 
will only weakly amplify the random field and hence lead to a weak 
change in polarized emission. 

We conclude that the roughly constant average polarized emis- 
sion across the arms and interam region cannot be easily explained 
with simple models of shock compression of magnetic field, if one 
simultaneously considers the weak contrasts that are observed in 
the total emission. The arm-interarm contrast in the radio emission 
probably results from a complex interplay of compression and de- 
compression of the dominant random field, occurring as the ISM 
undergoes phase changes on its passage through the arms. In addi- 
tion the thickness of the compressed regions compared to the width 
of the beam likely plays a role. 



8 SUMMARY AND CONCLUSIONS 

• Polarized emission (PI) is present throughout most of the disc 
of M51. In some regions the strongest PI coincides with the loca- 
tion of the strong spiral arms as seen in CO emission. In other re- 
gions PI is concentrated in the interarm region, forming structures 
up to ^ 5 kpc in size, reminiscent of the 10 kpc-scale magnetic 
arms observed in NGC 6946 (Beck & Hoemes,1996j ). The origin 
of these magnetic arms is still unknown. 

• The observed polarization angles trace spiral patterns with 
pitch angles similar to, but not always the same as, the gaseous 
spiral arms. The apparently ordered large-scale magnetic field re- 
sponsible for the well aligned polarization angles does not pro- 
duce a systematic pattern in Faraday rotation, leading us to con- 
clude that a large fraction of the polarized emission is caused by 
anisotropic small-scale magnetic fields (where small-scale refers 
to sizes smaller than the beam, typically 300-600 pc): anisotropic 
random fields, whose anisotropy is caused by a large scale pro- 
cess (for example, compression and/or shear) and so is aligned over 
large distances, can produce well ordered polarization angles but a 
random Faraday rotation distribution. 

• Faraday depolarization, caused by Faraday dispersion due to 
turbulent magnetic fields, leads to strong depolarization of the 
A20 cm polarized emission from the disc. From the observed fluc- 
tuations in Faraday depolarization we were able to estimate a typi- 
cal diameter of a turbulent cell as ~ 50 pc. 

• Fourier filtering followed by averaging in sectors is necessary 
to reveal the contribution of the regular (or mean) magnetic field to 
the observed polarization signal. This allowed us to fit a model of 
the 3D regular magnetic field of M5 1 to the observations of polar- 
ization angles at A A3, 6, 18, 20 cm. Due to the strong depolariza- 
tion at A A 18, 20 cm we were able to identify two different regular 
magnetic field patterns. In the thermal disc the regular field can be 
described as a combination of m = + 2 azimuthal modes, with 
the m = mode being the strongest: this combination can be the 
result of the strong two-armed spiral pattern modifying a dynamo 
generated m = mode (the easiest to excite according to mean- 
field galactic dynamo theory). The pitch angle of the m = mode 
is similar at all radii. In the halo the observed polarization angles 
at AA18, 20 cm, whose emission from the thermal disc is heavily 
depolarized, reveal a m = 1 azimuthal geometry for the regular 
magnetic field. The origin of this halo field is unclear. 

• The arm-interarm contrast in gas density and radio emission 
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was compared to a model where a regular and (isotropic) random 
magnetic field is compressed by shocks along the spiral arms. We 
found that where a strong arm-interarm contrast in PI is observed, 
in the inner arms r < 1.6 kpc, amplification of the random mag- 
netic field by compression can successfully explain the data, pro- 
vided that the regular magnetic field is not significantly increased. 
This constraint is similar to that obtained for two barred galaxies in 
|Bec k et al. (2005 ), where it was proposed that the regular magnetic 
field de-couples from molecular clouds as they rotate and collapse. 
We were unable to explain the average arm-interarm contrast in 
total and polarized radio emission, typical for much of the galaxy 
at r > 2 kpc, by a model involving shock compression of mag- 
netic fields. Even when the regular magnetic field is not compressed 
some increase in PI in the arms is expected from compression of 
the random field, whereas the average arm-interarm contrast in PI 
is about one: this problem could be resolved if the random field is 
isotropic in the arms and becomes anisotropic due to decompres- 
sion as it enters the interam. Alternatively, the compressed region 
of magnetic field may be sufficiently narrow (with a width of about 
20 pc), as might occur if the molecular clouds are de-coupled from 
the synchrotron emitting gas, to reduce the arm-interarm contrast 
by the required degree. 
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APPENDIX A: PARAMETERS OF THE FITTED 
REGULAR MAGNETIC FIELDS 

In Table |Al| we give the parameters of the fitted regular magnetic 
field models discussed in Section [6] Although a component of the 
regular field perpendicular to the disc plane {Rz in Eq. ( |A1| )) is 
allowed in the model and we searched for fits using this component, 
a vertical field was not required to obtain a good fit in either ring. 
The greater of the standard deviation and the noise within a sector 
was taken as the error in polarization angle. 
The regular magnetic field is modelled as 

Rr — Ro sinpo + sinpi cos((/) — /3i) 

+ R2 sinp2 cos(2(/) — /32), 
Rcf) = Ro cospo + ^1 cospi cos((/) — /3i) 

+ R2 C0SP2 cos(20 - /32), (Al) 

Rz = RzO + Rzl C0S((/) - Pzl) -\- Rz2 C0S(2(/) - l3z2), 

Rhr = Rm sinp/^o + Rhi sinp/^i cos((/) - hi), 
Rh^ = Rho cospho + Rhi cosphi cos((/) - /3hi), 

where Ri is the amplitude of the i'th mode in units of rad m~^, pi 
is its pitch angle, /3i (i ^ 1) determines the azimuth where the cor- 
responding non-axisymmetric mode is maximum and the subscript 
h denotes the components of the halo field. The magnetic field in 
each mode of this model is approximated by a logarithmic spiral. 
Pi — const, within a given ring. However, the superposition of such 



modes with different pitch angles leads to deviations from a loga- 
rithmic spiral. For further details of the method, see |Berkhuijsen et| 
|ari< [T997l l and |Fletcher etari ( |2004l ». 

The foreground Faraday rotation due to the magnetic field of 
the Milky Way R^g was also included in the fitting; we expect Rig 
to be the same in all rings and this provides a useful consistency 
check to the results of independent fits to the four rings. A log- 
arithmic spiral pattern has the same pitch angle in all rings. S is 
the residual of the fit and Sx the residual at a given wavelength. 
The appropriate value, at the 95% confidence level, is shown 
for the number of fit parameters and data points. A fit is statisti- 
cally acceptable if 5* ^ ^nd the Fisher criteria, that tests if S is 
unduly influenced by a good fit to a single wavelength, is satisfied. 
The values vary from ring to ring (even when the same number 
of fit parameters are used) as some sectors are excluded from the 
model, either because the noise in the sector exceeds the standard 
deviation of the measurements or because the sector represents an 
outlier from the global pattern. 

In Figs. |A1| |A2| |A3| and |A4| we show the observed sector- 
averaged polarization angles and the fitted model for each ring. 
The fit quality is excellent at A3 cm and A6 cm for all rings, but in 
the inner two rings sharp discontinuities in the A20 cm polarization 
angles around ^ 120° cannot be accommodated by the model. 
Since the parameters of the fitted halo field are largely determined 
by the data at the longer wavelengths we are therefore less satisfied 
with the parameters of the fitted halo field in these rings than with 
those of the disc field, as we discuss in Sect. [6] 



APPENDIX B: COMPRESSION OF A PARTIALLY 
ORDERED MAGNETIC FIELD IN A SPIRAL ARM 

We introduce a Cartesian frame in the sky plane {x,y,z) centered 
at the galaxy center with the x-axis directed towards the western 
end of the major axis and the ^-axis, in the northern direction; the 
z-axis is the directed towards the observer (and in the general direc- 
tion of the galaxy's north pole). We also introduce galaxy's Carte- 
sian frame {x,y, z) where the x- and x-axes coincide and the z- 
axis is also directed towards the galaxy's north pole. Magnetic field 
components in the two frames are related by (Berkhuijsen et al.| 
|T997l ): 

Bx — Bx, 

By — By COS i -\- Bz s'mi, (Bl) 
B2 = — By sini -\- Bz cos i, 

where i is the galaxy's inclination angle, and we include Bz for the 
future convenience (we recall that we can adopt 5^ = in M51). 
The galaxy's cylindrical frame r, 0, 2; then has the azimuthal angle 
measured counterclockwise (along the galaxy's rotation) from the 
X-axis. And finally, we introduce a local Cartesian frame of the 
spiral shock {x ,y' , z'), with the x^-axis perpendicular to the shock 
and directed from the interarm region into the spiral arm, the y'- 
axis parallel to the shock, so that the z^-axis complements them to 
a right-handed triad i.e., is directed towards the north pole of the 
galaxy. Then the angle between the x- and x-axes is 

= ct>-p.. (B2) 

It is convenient to specify the upstream large-scale magnetic 
field in the galaxy's cylindrical frame B = B{smpB, cosps, 0), 
where we neglect the vertical field component (see Sect. [6]). Then 
the unit normal to the shock, in the galaxy's frame is given by 
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Table Al. Parameters of fitted model for M51, with notation as in Eg. |A1| The index h refers to the "halo" field. The residuals at the wavelengths 
AA3, 6, 18, 20 cm are given in brackets. 







2.4 < r < 3.6 kpc 


3.6 < r < 4.8 kpc 


4.8 < r < 6.0 kpc 


6.0 < r < 7.2 kpc 




rad m ^ 


4 ±2 


5 ±4 


2 ±5 


4±i 


Ro 
Po 

R2 
P2 

^2 


rad m"^ 
deg 

rad m"^ 

deg 

deg 


-46 ±3 
-20 ±1 
-33 ±2 

-12 ±2 
-8 ±5 


^' -20 

-24 ±4 
-25 tl 
16 til 
-6 ±9 


' ^ -27 

-22 ±4 
-401- 

81^2 
-14 ±10 


-76 ±2 
-18±i 
-44 ±2 

3 ±2 

-25 ±2 


RhO 
PhO 

Rhi 
Phi 
Phi 


rad rQ— ^ 
deg 

rad rQ— ^ 

deg 

deg 


23 ±6 
-43 ±13 
76 ±11 
-45 ±5 
44 ±5 


77 ±26 

-49 ±12 
30 ±14 


57 ±26 

-50lf3 

o +23 
-29 


43 ±14 
-90 ±3 
-16 ±3 


Sx 
S 


(AA3,6,18,20cm) 


(14,21,21,21) 

77 

79 


(9,9,20,28) 

66 

83 


(10,18,11,16) 

55 

83 


(10,26,21,21) 

78 

79 



One data point excluded in the ring 2.4 < r < 3.6 kpc at: A20 cm, cf) = 140° 
Three data points excluded in the ring 6.0 < r < 7.2 kpc at: A3 cm, (/> = 20°; A18 cm, cf) = 220° 



A20cm, (/) = 180°. 



n = X ^ = (cos 0, sin 0, 0), the tangent vector is given by t = y' = 
( — sin ^, cos ^,0), and the regular magnetic field components nor- 
mal and tangent to the shock then follow, respectively, as 



where (6^,) 



(n) 



B^f = B sin(pb - Pa), Byf = B cos(p6 



Pa 



It is now easy to see that the compressed large-scale magnetic 
field in the spiral arm has the components 



where Cn = n^"^^ /n^^^ is the gas density compression ratio. Now 
we can transform both the original and compressed fields to the 
galaxy's frame by rotating it by the angle —0, Bx = B^' cos — 
Byf sin ^, By = Bx' sin + By/ cos ^, Bz = Bz' , and then to the 
sky frame using Eq. ( |B1| ), where the x and y component contribute 
to the magnetic field in the plane of the sky and is directed along 
the line of sight: 

B^ — (Bx,By), B\\ = Bs. 

Similar relations can be written for the random magnetic field 
6, but now we cannot neglect the ^-component of the random mag- 
netic field in the galaxy frame. The compressed random field in the 
shock frame is given by 

which can be transformed to the galaxy's frame and then to the sky 
frame to obtain 



bx — bx = bx' cos - Cnby' sin 6, 
by — 6^' sin^cosi + en^y' cos^cosi + Cn^z' sini, 
bz — —bx' ^1^0 sin i — tnby' cos siYii -\- Cnbz' cos i. 
Averaging using these equations then yields 





1 2(u) 


[1 + (4 


- 1) sin^ 0] , 




1 2(u) 


[1 + (4 


— 1)( — 1 sin^ 9 cos^ i)] , 


mr' = 


1 2(u) 


[1 + (4 


— 1)( — 1 sin^ 9 sin^ i)] , 



= icr^^"^ (isotropy of the 
upstream random magnetic field), and {bx'byf)'^^^ — {bx'bz')'^^^ — 
{by'bz')^^^ — (statistical independence of the upstream field 
components). 

Now the arm-interarm contrasts in various observables can be 
estimated as follows. For the total synchrotron intensity, 

^P2 (d) 

ei 



7 



where Lj is the pathlength through the synchrotron layer and 
the cosmic ray number density. The polarized emissivity in a par- 
tially ordered, anisotropic random magnetic field and uniform cos- 
mic ray distribution can be calculated using Eq. (20) of |Sokoloff et| 
[aLl pok] ) as 

PI oc {Bl - + (bl) - {bl)f + ABlBl 

Applying this formula to the downstream magnetic field, we obtain 
the contrast in polarized intensity epj by dividing it by PI^^^ oc 
B\'^^\ and also multiplying by the ratio of pathlengths Lj^^/Lj^-* 
and cosmic ray densities 

The Faraday rotation measure can be calculated as RM ^ 
O.SlneB^LRM, where Lrm is an appropriate pathlength, and its 
standard deviation is obtained using Eq. ([5]), and then the arm- 
interarm contrast is calculated straightforwardly. 
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